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Abstract 

The receive field of MRI imparts an image contrast which is spatially fixed rela- 
' r-| '- tive to the receive coil. If motion correction is used to correct subject motion oc- 

Q-i curring during an EPI time series then the receiver contrast will effectively move 

relative to the subject and produce temporal modulations in the image amplitude. 
This effect, which we will call the RFC-MoCo effect, may have consequences in 
the analysis and interpretation of fMRI results. There are many potential causes 
of motion-related noise and systematic error in EPI time series and isolating the 
" ^ . RFC-MoCo effect would be difficult. Therefore, we have undertaken a simu- 

^! lation of this effect to better understand its severity. The simulations examine 

P_J this effect for a receive-only single-channel 16-leg birdcage coil and a receive- 

only 12-channel phased array. In particular we study: (1) The effect size; (2) 
Its consequences to the temporal correlations between signals arising at different 
spatial locations (spatial-temporal correlations) as is often calculated in resting 
state fMRI analyses; and (3) Its impact on the temporal signal-to-noise ratio of 
an EPI time series. We find that signal changes arising from the RFC-MoCo 
effect are likely to compete with BOLD (blood-oxygen-level-dependent) signal 
changes in the presence of significant motion, even under the assumption of per- 
fect motion correction. Consequently, we find that the RFC-MoCo effect may 
lead to spurious temporal correlations across the image space, and that temporal 
SNR may be degraded with increasing motion. 
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1. Introduction 



Subject moti on can be a significan t source of noise in BOLD (blood-oxygen- 
level-dependent (|Ogawa et all 1 19901) ) functional magnetic resonance imaging 
(fMRI), a neuroscience research tool employing time series of gradient echo 
images, most commonly using the echo planar imaging (EPI) pu lse sequence . 
This problem has been recognized since the inception of fMRI (IHajnal et all 
1994h . The most obvious concern is that motion of the brain, with its contrast 
fixed relative to the head, will produce temporal variations at any given point in 
image space that can mask the BOLD or be erroneously attributed to the BOLD 
effect. But there can also be contrast within brain images that is fixed relative to 
the scanner. Such contrast can potentially be exascerbated by motion correction 
since following correction such contrast would then move relative to the brain 
and produce its own temporal variation. 

With the proliferation of MRI systems employing multichannel receiver ar- 
rays (phased arrays) arises greater potential sensitivity to motion-related error 
and noise in EPI time series data due to scanner- fixed receiver contrast. Although 
multichannel receiver arrays hold great promise for improved tSNR (temporal 
signal-to-noise ratio) it presently appears th at motion-related noise may be lim- 
iting their usefulness (IHartwig et all 1201 1|) . The effects of motion can be espe- 
cially problematic when acceler ated imaging (IGriswold et al.Ll2002l : |Pruessmann et al 
1999|) is used in EPI time series (|Sheltraw et all 1201 2|) . However, even in the ab- 
sence of acceleration there are many potential motion-related effects that can 
prevent multichannel receiver arrays from achieving their motion-free tSNR po- 
tential. One such effect, arising from scanner-fixed contrast, which we will call 
the RFC-MoCo effect (Receive Field Contrast Motion Correction effect), may 
occur when motion correction is applied to EPI time series data possessing sig- 
nificant receive field contrast. 

There are many potential causes of motion-related noise in EPI time series. 
Here are some of the acknowledged important sources of noise and systematic 
error in non-accelerated EPI: 



1. Poor motion correction: Most motion correction algorithms used in fM RI 
(Eox[ll996l:lFriston et all[l995l:ISmith etalll200ll : IWoods et all[l99ll) as- 
sume that a volume of image slices moves as a rigid body. Since a volume 
of data is usually collected over a time interval of 1 to 2 seconds there is 
clearly enough time for the head to move during such an interval, thereby 
violating the rigid volume assumption. This may lead to bad estimates 
of motion parameters and inadequate correction of motion-related contrast 
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fluctuations. Furthermore the issue is often complicated by the use of addi- 
tional temporal interpolation (lRochel . 1201 luBannister et al.Ll2007h as well. 

2. Spin history: The steady state of the pre-excitation magnetization, estab- 
lished after a few TR (repetition time) periods in the absence of motion, is 
perturbed by motion o rthogonal to the imaging planes of 2D multislice EPI 
(IMuresan et al.Ll2005h . potentially causing increased noise and systematic 
efl'ects upon spatial-temporal correlations. 

3. Main field inhomogeneity: The presence of an object within the strong 
main magnetic field Bq of MRI leads to increased inhomogeneity of the 
field. The inhomogeneity r esults in image distor t ions w hich increase in 
severity with field strength (IJezzard and BalabanL Il995|) . Shim coils are 
used to try to restore field homogeneity but shim coils are fundamentally 
limited with respect to the fields they can compensate. When the object 
moves the inhomogeneity varies with respect to time which may poten- 
tially increase noise and systematic error in the time series images while 
decreasing image-space signal strength. 

4. Receive field distortio n: A receiver coil loade d with a head will distort 
the coil's receive field (IGiovannetti et al.L 12008). This distorti on should be 



more pronounced at high main field strengths (|Wiesinger et al.l . 120061) . As 



the head moves this distortion will vary in time thereby adding temporal 
fluctuations in an EPI time series. 
5. RFC-MoCo effect: The receive field, when not spatially hom ogeneous, 
will i mpose an image contrast that is fixed relative to the scanner (|Hartwig et al. . 



201l|) . When an assumed perfect motion correction is applied to a time 



series of such images the receive field effectively moves relative to the 
imaged object. This temporally varying image contrast may potentially 
produce tSNR degradation and systematic effects upon spatial-temporal 
correlations. 

6. TFC-MoCo effect: The transmit field, when not spatially homogeneous, 
will also impose an image contrast that is fixed relative to the scanner and 
will produce temporally varying image contrast in a motion-corrected time 
series of images in a manner analogous to the RFC-MoCo effect. Most 
commercial scanners used in fMRI research today make use of a large 
diameter birdcage body coil to produce the transmit field. The large diam- 
eter, as compared to a head coil, increases the transmit field homogeneity 
and therefore decreases the potential TFC-MoCo effect. 

The RFC-MoCo effect is the focus of this paper. In particular, we seek to es- 
timate the effect size of tSNR degradation and any spurious spatial-temporal 
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correlations induced by the receive field contrast. This is an important subject 
because there is a small handful of reports of little benefit arising from use of a 



32-channel coil for fMRI, compared to a 12-channel coil (IHartwig et all |201 1 



Kaza et all I2009L l2011l : iLi et all l2009|) . and we have received anecdotal reports 



of poorer performance of the 32-channel coil for fMRI in our own center. Thus, 
we were motivated to investigate the RFC-MoCo effect because of a suspicion 
that it may produce a limiting motion sensitivity when compared to the same 
protocol acquired with a smaller phased a rray coil. Furthermore, recent reports 



(IVan Dijk et al.L 1201 2l : iPower et al.l . 120111) have suggested that a significant (but 



unknown) fraction of correlations determined in resting-state functional connec- 
tivity (RS-FC) studies are driven by head motion. We hypothesized that the 
RFC-MoCo effect could be a cause of such errors, given the use of phased array 
coils for the majority of RS-FC fMRI studies to date. 

The use of a prescan normalization as a means for correctin g the RFC-MoCo 
effect has been stu died by Hartwig et al. (|Hartwig et al.L 1201 1|) . Also Kaza et al. 



(IKaza et al.L |2009|) used prescan-normalized and unnormalized data in a com- 



parison of fMRI efficacy using 12-channel and 32-channel phased array coils. 
However, experimentally isolating the RFC-MoCo effect from the other potential 
causes of motion-related effects in time-series EPI would be very difficult if not 
impossible. Therefore in this paper we undertake simulations of the RFC-MoCo 
effect, which by design incorporates an assumption of perfect motion correction, 
that will allow us to elucidate the effect's impact upon EPI time series data. 

When choosing a receive coil for fMRI it is important to assure that the sub- 
ject within the coil will have visual access to stimuli or cues presented on a screen 
placed outside the coil. Commercially available 12-channel arrays give good vi- 
sual access for fMRI applications and are widely used. A birdcage coil with 
sixteen legs will give comparable visual access to that of the 12-channel array 
consisting of the usual overlapping coil elements. This specification motivates 
our choice of coil geometries in our simulations. In this paper we carry out sim- 
ulations of the RFC-MoCo effect for: (1) A 12-channel head coil receive array 
consisting of 12 independent coil elements with a cylindrical geometry similar 
to those commonly used in modern scanners; and (2) A 16-leg receive-only bird- 
cage head coil. In particular we will simulate: (1) The RFC-MoCo effect size; 

(2) A spatial map of temporal correlations due to the RFC-MoCo effect with re- 
spect to a seed point (as is often calculated in resting state fMRI analyses); and 

(3) The tSNR of an EPI time series when the noise of the RFC-MoCo effect is 
present in addition to the usual Gaussian noise of the complex-valued images. 
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2. Theory 



The use of multichannel receiver arrays is commonplace in magnetic reso- 
nance imaging and fMRI in particular. Whether the image data is to be used 
in accelerated GRAPPA-like parallel imaging or non- accelerated imaging, one 
must combine the image data from the individual coil elements into a single com- 
posite image. Of the various methods that could be used to generate a composite 
image from the image s generated by each coil, the sum-of-squares (SOS) method 



(IRoemer et al.Lll990h is ubiquitous. The SOS image reconstruction method cre- 
ates a composite estimated image p(r) from the complex- valued images p,„(r) 
associated with each of the M receive coils of the array by the following opera- 
tion: 



M 



m= 1 



1/2 



(1) 



Note that when M = 1 (a single channel) the result is the usual magnitude image 
associated with a birdcage coil. Equation ([T]) assumes that we have a represen- 
tation of the image p(r) in the continuum image- space rather than a discrete 
image-space. In practice the image will always be in a discrete image-space 
which is related to the underlying continuum image-space object through a con- 
volution with a point-spread function determined by the sampling of fc-space. 
The discretization is expected to influence the RFC-MoCo efi'ect upon tSNR, as 
well as methods to correct for the efi'ect, but for the purpose of establishing the 
size and importance of the efi'ect the use of the continuum image-space represen- 
tation should be sufficient and is expected to give a best-case estimate of the size 
and importance of the efii'ect in EPI time series used in fMRI. 

In the absence of noise (electronic or body noise) and motion Equation ([T]) 
can be written as 



P(r) 



M 



S \p{r)c 



m=l 



1/2 



M 

s 

m=l 



1/2 



(2) 



where p(r) is the complex valued true image and c„,(r) is the receive field of 
the m''' coil element of the array. Therefore the SOS method results in an image 
contrast which is dependent upon the array geometry and is spatially-fixed with 
respect to the scanner. 

The imaged object p may undergo rigid body motion (rotation and transla- 
tion) with respect to the scanner frame of reference. We denote this motion by a 
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time-dependent affine transformation operator J5^(?„) where t„ = nAt is the time 
at which the n''^ volume of image data is acquired. We will assume that the object 
does not move by a significant amount on a time scale less than At. This restric- 
tion imposes another best-case scenario on the effects of motion: Real motion is 
likely to occur within a TR period and will produce more complex effects, in par- 
ticular spin history effects that will interact with the RFC-MoCo effect. Thus, for 
simplicity, we consider the RFC-MoCo effect in isolation and note that the real 
effects on data are likely to be worse than presented here. The SOS estimated 
image of the object at time t,, is then given by 



P(r, tn) 



M 



2 \c,n{r)Jl{t„)p{r) + ?7„ 



1/2 



(3) 



where ?/„, is the uniform noise image for the m''' coil. 

The usual, but idealized, motion correction is performed by applying the 
inverse of the affine transformation Jl{tn) to the SOS image p(r, tn) to obtain the 
motion-corrected image pc(r, t„) given by 



Pc{r,tn) 



^-\tn)p{r,tn) 



M 



2 \Cn,{r)j{{tn)p{r) + T]„ 



m=l 



M 



2 |p(r)^ \tn)c„,{r) + r]„ 



m— 1 



1/2 



1/2 



(4) 
(5) 

(6) 



Therefore the effect of combined motion and motion correction is mathemati- 
cally equivalent to moving the receive coil elements relative to the imaged object. 
As explained in the Methods section this equivalence will be used to generate the 
time varying contrast of our simulated RFC-MoCo effect. H 

To calculate the RFC-MoCo effect size, temporal correlation map, and tSNR 
we respectively calculate the percent difference map Dp{r), correlation value 
mapxi^P^ ^p') tSNR map tSNR(r) defined according to: 



If necessary the action of the affine transformation operator upon an image can be exphcitly 
written as ^p{v) = p(Rv — To) where R, a matrix representing rotation about the object's center 
of mass, and To, a vector representing translation of the center of mass, are functions of f„. 

^If need be, the inverse affine transformation can be written more exphcitly as J?[^'p(r) = 
p(R-'(r-r„)). 
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tSNR(r) 



E[{pcirp,tn) - Hr„){pc{rp>,tn) - Hr.)] 



E[p,.(r,f„)] 



VE[|A-(r,r„)P]-|E[^c(r,r„) 



(7) 



(8) 



(9) 



where E denotes the expectation value with respect to time f„ while and cTr are 
respectively the mean and standard deviations of the simulated (r, t„) . 

When calculating the percent difference and the correlation map for our sim- 
ulations we will neglect the noise term 7/,„. This results in the following quantities 
which are independent of the imaged object p(r): 



E[{C{rp,tn) -yUr,J(C(v,?„) -lUr,,,)] 



(10) 



(11) 



where the net receive field contrast C(r, tn) is given by 



M 



1/2 



C{r,tn) = 2] |^"n^n)Cm(r) 



(12) 



vm=l 



When calculating the simulated tSNR we will assume that the image is uniform, 
p(r) = 1, so that we may isolate the RFC-MoCo effect from brain contrast 
effects. In this case we may then write 



tSNR(r) 



M 



1/2- 



\m— 1 



(13) 



M 



;n=l 



\ 



M 



J] \ JI ^{t„)c,„{r) + r]„ 



m — i 



When the motion of the object is a small translation, that is ro{t) = 6r{t) (see 
Appendix C[) we can approximate Equations (flOl) and (fTTI) by 




[6r{t2)-Sr{h)]-VC,o.{r,t2) 



X 100 



(14) 



E[{6r{tn) ■ VC,„,(r^,0))((5r(?„) • VC,„,(r^s 0))] 



(15) 



^E[{6r{tn) ■ VC,,,(rp,0))2] ^E[{6r{t„) ■ VC,,,(v,0))2] ' 



where in deriving Equation (fTSl) we have assumed that E[6r{tn)] = as will be 
the case in our simulations. 

Although Equations (fT?l) and (fT5l) will not be used in our simulations of the 
RFC-MoCo effect these equations can be helpful with respect to understanding 
the effect. For example Equation (fT4l) shows that the percent difference should 
be directly proportional to the translational displacement and the gradient of the 
SOS contrast field. Also from Equation (fTSl) it should be clear that if the small 
translations are in only one direction then the correlation map takes the value 1 
or —1 depending upon the sign of VCsos{^p) and VCsosi^p') the seed and map 
points respectively. When the motion in each dimension is temporally different 
then —1 ^ xi^p^'^'p') ^ 1 and the correlation map can have a more interesting 
structure, as will be seen in the results section of this paper. 

3. Methods 

In this section we: (1) Describe the calculation of the simulated temporally 
varying receive field contrast for the 16-leg birdcage and 12-channel array head 
coils; (2) Describe the calculation of the metrics of the RFC-MoCo effect - per- 
cent difference, temporal correlation and tSNR maps - in a representative trans- 
verse (orthogonal to the long axis of the head coil) plane at z = 90 mm. Each 
coil has a radius of 130 mm at the plane of interest and the metrics are calculated 
over alio mm radius region centered in this plane. This is equivalent to assum- 
ing a simulated spherical phantom having a diameter of 220 mm - comparable to 
a typical image field-of-view for an adult human head - which is centered within 
the head coil. A uniform signal intensity from all points in space was assumed, 
thereby allowing us to investigate the RFC-MoCo effect in the absence of image 
contrast, transmission field heterogeneity and magnetic susceptibility gradient 
effects. 
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3.1. Simulating the Temporally Varying RFC 

To calculate the receive fields of the head coils we construct the coil or coil 



elements as a composite of line segments of current. As shown in Appendix A 
the receive field due to a single line segment with specified endpoints and current 
amplitude can be written in an exact mathematical form by means of a Biot- 
Savart Law integration. Then to translate the receive field at time f„, thereby 
creating a time series of simulated RFC-MoCo data, we simply translate each 
endpoint of the line segments comprising the head coil or coil elements. All 
calculations were performed using C++ code and since the mathematical form 
of the receive field is known exactly then all calculations will be exact to within 
machine precision. Note that for both head coils we will neglect the efi'ects 
of coil loading and mutual inductance between current-carrying line segments. 
Although the error due to these approximations is significant in the precise design 
and operation of receive coils it should not change the results given in this paper 
significantly. 

Some of the simulations will be performed using a time series of realistic 
translational head motion data. Figure[I]shows this realistic center-of-mass trans- 
lation as a function of time. This motion data was obtained from the output of 
the FSL (FMRIB Software Library, Universi t y of O xford, UK.) motion correc- 



tion algorithm (MCFLIRT) (|Jenkinson et al.L |2002[) for a human subject in an 



actual fMRI experiment and is representative of motion data obtained on normal 
adult subjects at 3 T. The temporal mean of the motion was subtracted from the 
data in accordance with the assumptions of our motion model. To test the effects 
of variable overall motion the motion data is variably scaled to yield specific 
amounts of root-mean- square motion. 

16-Leg Birdcage Coil 

Figure [2] depicts a birdcage coil of the type used in this paper. In mode 1, the 



most spatially homogeneous mode (|JinL 11989^ the time-independent part of the 
mesh current /„ is given by: 

/„ = Ccos(2;m/A^) (16) 

where A'^ is the number of legs (struts) in the coil and C is a constant which will 
not be of consequence in this paper. The current J„ in the n"' leg is given by 

The endpoints of line segments comprising the legs and endrings of the bird- 
cage coil are arranged symmetrically on a cylindrical surface of 130 mm radius. 
Figure [3] shows the receive field for such a birdcage coil. The field is plotted over 
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Figure 1 : The time series of translational displacements used in the tSNR and correlation value 
calculations. 




Figure 2: The simulations are done with a 16-leg birdcage coil but here we depict a 12-leg coil to 
avoid clutter in the figure and clearly convey the general coil geometry. The coil is depicted with 
mesh currents /„ (« = 1 , . . . , A^) where is the number of legs. Note that the lower end ring lies 
in the z = mm plane. 
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an axially- sliced region of a 1 10 mm radius which is centered transversely within 
the coil. The important thing to note about Figure[3]is that the birdcage coil gives 
a receive field which varies by approximately 30% over the region shown. 




100 



100 



Figure 3; Receive field magnitude for 16-leg birdcage coil geometry (coil radius = 130 mm, 
coil length L = 186 mm). The field is plotted over a region of radius 110 mm centered in the 
transverse plane at z = 90 within the coil. The receive field magnitude is plotted as a unitless 
quantity since the homogeneity of the field is our primary concern in this paper 



12-Channel Array 

The simulated 12-channel array consists of M = 12 coil elements of the 
general form shown in Figure |4] in which the coil element is comprised of six 
segments each with a unit current in the clockwise direction from the perspec- 
tive of an observer external to the array. Figure [5] depicts segment 1 of each coil 
element, which are taken to lie in a single transverse plane (the x,y-plane). The 
endpoints of the six line segments comprising each coil element lie on a cylin- 
drical surface of 130 mm radius (see [Appendix B[ ) and closely approximate the 
geometry of commercially available head coils. The sum of the fields due to each 
of the line segments comprising a given coil element yields the receive field of 
that coil element. The final SOS receive field is then calculated from the receive 
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fields of the individual coil elements according to Equation Q. 



190 mm 80 mm 




Figure 4: The general shape of a single coil element for the 12-element array coil. The length of 
segments 2 and 6 are the same for all elements. Segments 1 lie in the z = transverse plane. 

Figure[6]shows the SOS receive field for the simulated 12-channel array. The 
field is plotted over an axially-sliced region of 1 10 mm radius which is centered 
transversely within the receive array. The important thing to note about Figure 
[6] is that the 12-channel array gives a receive field which varies by as much as 
400% over the region shown. 

3.2. Calculating RFC-MoCo Effect Metrics 

Percent Difference 

For both head coil geometries we estimate the magnitude of the RFC-MoCo 
effect size by calculating the percent difference map, given by Equation (flOl) . 
when the object is displaced translationally by 1 mm in the y-direction. For both 
head coils we also calculate the average and maximum of the percent difference 
map over the 1 10 mm radius region when the object is displaced by to 3.5 mm. 
For the percent difference calculations no noise is added to the simulated data. 

Temporal Correlation 

We examine potential effects upon temporal correlation maps due to the 
RFC-MoCo effect for translational head motion in the 12-channel array only. 
These simulations are performed by using the time series of realistic head motion 
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Figure 5: Cross-section of the simulated 12-channel receiver array showing length and orienta- 
tions for segment 1 of each coil element. The z-axis (parallel to Bq), of a right-handed coordinate 
system, points into the plane of the page. 




Figure 6: SOS receive field over an axially-sliced region of radius 1 10 mm which is centered in 
the transverse plane at z = 90 within the 12-element receive array. The receive field magnitude 
is plotted as a unitless quantity since the homogeneity of the field is our primary concern in this 
paper 
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variables (see Figure [B and calculating temporal correlation maps according to 
Equation (fTTI) . To calculate the temporal correlation maps we set the root-mean- 
square (rms) magnitude of the translational motion to 1 mm but it should be 
noted that for small translations we expect, from Equation (fT5l) . that the correla- 
tion maps will be independent of the magnitude of the motion. The correlation 
is between a seed point at ;c = —104 mm, y = mm and z = 90 mm and 
all points within the same axially-sliced (z = 90 mm) region of radius 110 mm 
which is centered within the array. For the temporal correlation calculations no 
noise is added to the simulated data because our principal concern is elucidating 
the general features one might expect to see in a temporal correlation map when 
the RFC-MoCo effect is of significant sizeH. 

tSNR 

For the 12-channel array only we investigate the spatial dependence of the 
tSNR maps for varying amounts of realistic translational head motion (see Figure 
[T]) and spatially uniform Gaussian electronic noise 77,,,. The noise is applied to the 
complex-valued time series image data according to Equation (fTSl) such that it 
yields a tSNR level, 80 at the center of the image, representative of the tSNR seen 
in a 3-minute 72* -weighted EPI time series at 3 T. TSNR maps were obtained 
for rms motion values of 0.0, 0.5, 1.0 and 2.5 mm. 

4. Results 

Percent Difference 

Figures |7] and [8] show the percent difference in the receive field contrast due 
to a 1 mm translation in the ^'-direction for the 16-leg birdcage and 12-channel 
array head coils, respectively. From these figures it should be clear that this 
particular birdcage coil geometry yields significant reduction in the RFC-MoCo 
effect as compared to the 12-channel array. It is also clear that the RFC-MoCo 
effect for the 12-channel array is of similar magnitude if not greater than the 
percent difference expected from the BOLD effect at most points in this axial 
slice. Hence the effect should be an important systematic error in fMRI analysis 
if an array coil is used as the receiver. 

Figure H] shows a plot of the average percent difference and the maximum 
percent difference over the 110 mm radius region of an axial slice at z = 90 mm 
versus the displacement in the y-direction. The trend is approximately linear, as 
would be expected for small displacements (see Equation (IC.4I) ). 



The significance is established by the percent difference maps 
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Figure 7: Percent difference in receive field contrast due to 1 mm motion in the y-direction for 
tfie 16-leg birdcage receive coil geometry (radius =130 mm). The percent difference is plotted 
over an axially-sliced region of radius =110 mm which is centered within the array. 




Figure 8: Percent difference in receive field contrast due to 1 mm motion in the y-direction for 
the 12-channel receiver array geometry (radius = 130 mm). The percent difference is plotted 
over an axially-sliced region of radius =110 mm which is centered within the array. 
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Figure 9: Comparison of percent difference between 16-leg birdcage and 12-channel array re- 
ceiver head coils. Average and maximum percent difference over an axially-sliced region of 
radius = 110 mm which is centered within each coil are plotted versus displacement in the y- 
direction. 



Temporal Correlation 

Figure [10] is a surface plot showing the temporal correlation value as a func- 
tion of location within an axial slice. The correlation is between a seed point at 
X = — 104 mm, y = mm and z = 90 mm and all points within the same axially- 
sliced (z = 90 mm) region of radius 110 mm which is centered within the array. 
For reference the value ofx at the seed point (1.0 as expected) is indicated by 
a red dot on the 3D surface. Notice the negative as well as positive correlation 
values. 

tSNR 

Figures [TT] through [T?] show the tSNR surface plots for the cases of no mo- 
tion, 0.5 mm, 1.0 mm and 2.5 mm rms motion for the 12-channel receive array. 
Color and viewpoint is different for each surface plot to aid in the visualization 
of the surface. Please refer to the axes for quantitative information. At 0.5 mm 
rms motion (Figure [T2l) much of the tSNR benefit from the multiple channels is 
eliminated through the RFC-MoCo effect. At 1 .0 mm rms motion (Figure 
the flattening of the tSNR map is effectively complete. At 2.5 mm rms motion 
(Figure [T4l) . which approximately corresponds to a pixel shift for the chosen 
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Figure 10: A surface plot showing the temporal correlation value xi^p, r^/ ) as a function position 
within an axial slice. The correlation is between a seed point Tpi at x = —104 mm, y = mm 
and z = 90 mm and all points Tp within the same axially-sliced (z = 90 mm) region of radius 
= 110 mm which is centered within the array. For reference the value of x at the seed point is 
indicated by a red dot. The motion was purely translational (1 mm rms) and taken from the ouput 
of the FSL (FMRlB's Software Library) MCFLIRT motion correction code (see Figure[T]i. 

electronic noise tSNR, the tSNR in regions nearest the coils is much reduced 
compared to that at the center of the array. Degradation of tSNR would be ex- 
pected for all brain regions. 

5. Discussion 

We set out to assess how receive field contrast could possibly limit the fidelity 
of motion-corrected time series MRI measurements acquired in the presence of 
significant subject motion. Simulations were performed so that the specific in- 
teraction of the receive field contrast and motion could be assessed without the 
complicating factors, such as head-fixed image contrast and imperfect motion 
correction, that would likely arise in experimental data. 

We have shown in this work that: (1) The RFC-MoCo efl'ect should be ex- 
pected to compete with the BOLD efi'ect for a typical 12-channel cylindrical 
array, (2) Interesting temporal correlations occur as a result of the RFC-MoCo 
effect that will compete with the temporal correlations due to the BOLD ef- 
fect and (3) Moderate amounts of motion lead to serious tSNR degradation in 
regions of the 12-channel array receive field where, from consideration of the 
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tSNR 




Figure 1 1: A surface plot showing the tSNR for the case of no motion as a function of location 
within an axially-sliced (z = 90 mm) region of radius =110 mm which is centered within the 
12-element array. The tSNR is 80 at the center of this region. 



240 r 



tSNR 




Figure 12: A surface plot showing the tSNR for the case of 0.5 mm rms motion as a function 
of location within an axially-sliced (z = 90 mm) region of radius = 110 mm which is centered 
within the 12-element array. The motion was purely translational and taken from the ouput of 
the FSL (FMRIB's Software Library) MCFLIRT motion correction code. 
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tSNR 




Figure 13: A surface plot showing the tSNR for the case of 1.0 mm rms motion as a function 
of location within an axially-sliced {z = 90 mm) region of radius = 110 mm which is centered 
within the 12-element array. The motion was purely translational and taken from the ouput of 
the FSL (FMRIB's Software Library) MCFLIRT motion correction code. 



tSNR 




Figure 14: A surface plot showing the tSNR for the case of 2.5 mm (the size of a pixel typically 
associated with the tSNR at the center of the array) rms motion as a function of location within an 
axially-sliced (z = 90 mm) region of radius = 110 mm which is centered within the 12-element 
array. The motion was purely translational and taken from the ouput of the FSL (FMRIB's 
Software Library) MCFLIRT motion correction code. 
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signal magnitude alone, one might expect improved tSNR. 

Caution should be exercised in interpreting the results given here for several 
reasons: 



1 . The effects of mutual inductance and coil loading were both neglected. We 
do not expect the neglect of these two complicating effects to do anything 
but create greater local variability of the receive array SOS pattern and 
hence even greater RFC-MoCo effect. 

2. We have assumed the current-carrying conductors of the simulated head 
coils to be thin wires. The conductors of real head coils are thin tracings 
with a width (parallel to the cylidrical surface of the head coil) in the neigh- 
borhood of 3 mm. The width of the real conductors should have a spatially 
smoothing effect which would be expected to reduce the RFC-MoCo effect 
to some extent. 

3. We used a root-sum-of-squares reconstruction of the final image. While 
this method of reconstructing images from coil arrays is ubiquitous, other 
image combination methods are avail able such as the adaptive combina- 
tion of coil images dWalsh et al.L 120001) . and these methods should perform 
differently with respect to the specific results of RFC-MoCo. 

4. We have not simulated rotational motion of the head in this work and the 
presence of rotation should be expected to give even more interesting cor- 
relation maps. 

5. Some manufactures combine the receive fields from the 12 coil elements 
in various manners to produce arrays that may be effectively smaller (eg. 
the "CP" and "Dual" receive modes of Siemens 3T Trio scanners) with 
correspondingly more homogeneous receive fields. These coil combina- 
tion approaches should establish reduced RFC-MoCo effects compared to 
the full multi-element array operation, but the remaining RFC-MoCo ef- 
fect is likely to be larger than for a birdcage coil having similar physical 
dimensions. 



In a real fMRI experiment our assumption of perfect motion correction will 
not be satisfied and the degree to which the RFC-MoCo effect will influence the 
time series data will depend upon the accuracy of the applied motion correction 
algorithm. When imperfect motion correction is applied there will be a mixing 
of the various contrasts present in the image - those fixed relative to the head 
and those fixed relative to the scanner. Furthermore, the presence of receive 
field contrast, as with other contrasts that are fixed relative to the scanner, is 
expected to degrade motion correction performance through an underestimate of 
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the actual motion - an effect we term the "anchoring effect." We are investigating 
the magnitude of the RFC-anchoring effect in parallel work. (Other scanner-fixed 
contrast mechanisms will generate their own anchoring effects.) 

The question naturally arises as to the practical relevance of the results pre- 
sented here for the RFC-MoCo effect. In the first instance, the 12-channel array 
was simulated in a manner that may not be appropriate for some scanner vendors. 
For example, in the absence of accelerati on the default coil reception mode on 
a Siemens scanner would be "CP mode" (IReykowskil |2006|) which is designed 
to provide near optimum SNR in the center region of the image. Thus, we can- 
not yet give a clear ranking of the RFC-MoCo effect's contribution to spurious 
correlations due to all causes in resting-state fMRI studies. We will attempt to 
clarify this ranking in future simulation work. 



Data processing methods such as independent component analysis (IKundu et al.L 



2012h may be able to discriminate between artifactual and neurally-driven corre- 



lations in some circumstances, but it has not yet been established that the motion 
metrics output from the commonly used affine correction algorithms will capture 
sufficiently all t he correlations intro duced by the RFC-MoCo effect. Recent work 
by Power et al. ([Power et al .11201 11) suggests that some spurious correlations can 
be disambiguated by removing networks correlated with motion parameters, but 
further work is needed to establish the efficacy of this approach. 

In follow up simulations we will provide a similar analysis to that presented 
in this work for the CP, dual, triple modes of the 12-channel head array. We 
will also present a similar analysis for a 32-channel array which has smaller 
somewhat hexagonally shaped coil elements arranged on the surface of a roughly 
spherical shell. Compared to the results given in this work we would expect the 
RFC-MoCo effect to be larger for the 32-channel array and smaller for dual and 
CP modes of the 12-channel array. We will then be in a better position to make 
quantitative estimates of the severity of the RFC-MoCo issue for connectivity 
studies from resting state data, and to compare the relative performance of 32- 
channel and 12-channel array coils for time series EPI under motion limiting 
regimes. 
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{xo,yo^zo) 



X 



Figure A. 15: A line element specified by the endpoints (xo,3'o,Zo) and {x\,y\,z\) and carrying 
a unit current in the direction of the unit vector e/. The receive field is to be calculated at points 

Appendix A. Magnetic Field Due To An Arbitrarily Oriented Current Car- 



In this section we derive the receive field due to a single line element. This 
result will be used to calculate the receive field of a coil element that can be 
constructed from a set of such line elements. Figure IaTsI depicts one such line 
element. 

The magnetic field B due to the line element can be calculated from the 
following form of the Biot-Savart Law: 



rying Line Segment 




(A.l) 



The unit vectors 6/ and e,. are given by 




(A.2) 



where Ax = x\ 



xq etc and 




X + 



y - y' ^ z- z! 

r r 



z 



(A.3) 



r 



where 




(A.4) 
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and 



Since 



then 



r = \ {x- x'Y + {y- y'Y + {z- z'f 



X = Xq + Iti ■ X = xq + {l/L)Ax 
y =yo + lery = yo + {l/L)Ay 
z' = Zq + lei ■ z = zo + {l/L)Az 



Gr = r ^[x — Xo — {l/L)Ax] X 
+ r-'[y-yo-{l/L)Ay]y 
+ r-'[z-zo-{l/L)Az]z 



and 

r = 
Since 



6/ X 



X y z 



L L L 

x-x' y-y' z-z' 



= (Lr) 

+ (Lr) 

+ (Lr) 

= (Lr) 

+ (Lr) 

+ (Lr) 



-'[Ay{z-z')-Aziy-y')]x 
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-'[Axiy-y')-Ayix-x')]z 
-'[Ay{z zo lAz/L) Az{y yo lAy/L)] x 
[Az{x - Xq- lAx/L) - Ax{z - Zo - lAz/L)] y 
'^[Ax{y ~yo~ lAy/L) - Ay{x - Xq- lAx/L)] z 
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(A.9) 
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where 



= {x- xof -2{l/L)Ax{x- xo) + il/Lf{Axf 

+ {y- yof - 2{l/L)Ay{y yo) + {l/L)\Ayf 

+ {z-ZoT-2{l/L)Az{z-Zo) + {l/L)\Azf (A.19) 

= (jc - xof + {y- yof + {z- zof 

- 2{l/L)[{x ~ xo)Ax +{y- yo)Ay + (z - zo)Az] + l\ (A.20) 

Substituting 

n( \ Ay z_zo_lAz/L Az C'^y_yo_lAy/L,, 

BJx,y,z) = — — TTTT"' tt; tttt"' (A. 21) 
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where 



^ = (x - xo)' + (); - );o)' + (z - zo)' (A.24) 
= {2/L)[{x-xo)Ax+{y-yo)Ay+{z-Zo)Az]. (A.25) 



From the integral tables 
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We can then write 



Br = ^ 
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(A.28) 



(A.29) 



(A.30) 



(A.31) 



Note, as a consistency check, that Equations (IA.28I) . (IA.30I) and (IA.31I) reduce to 
the fi eld due to an in finite length wire with unit current oriented along any of the 
axes (lArfkenl . Il985h . The receive field for a single line segment wire with unit 
current directed from point {xo,yo,zo) to point {xi,yi,zi) is then 



B* = Bx — iBy 
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Appendix B. 12-Channel Array Geometry 

The first end point (^°„,,y°„) and the second end point {x\^,y\^ of m''' seg- 
ment of the n'^ coil element are: 

Top half of array 





= r 


y\:- 


= 




= r cos 01 


= r sin 01 
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yVi 


= — rsin(02 + (p\ — 62 — 


e.) 



where 0„ = 2 sin^' (/„/2r) is the angle subtended by a chord of length /„ giving 
the length of segment 1 of the n'^ coil element and 9n = 2 sin~^(J„/2r) is the an- 
gle subtended by a chord of length giving the chord length of the coil overlap 
and gaps. 

Bottom half of array 
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where 0, the angle subtended by segment 1 for coils 5 through 8, is determined by 
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the condition TT- 2 (04 + ^0-^0 = 40-3^2 or0 = ^4- (04 + ^0 -^i)/2 + 36*2/4. 
We choose r = 130, do = l,di = 20, J2 = 23, and Zi = 95, = U = /lo = Zg = 
77. For segment 4 we choose r = 50 with all the same angles as for segment 1. 
See Figure |4] for the z-coordinates of each endpoint. 

Appendix C. Small Rotations and Translations 

Let /(r) be a scalar function of the vector r and let r' = R^(r — r^). For 
infinitesimal rotations and translations we can write r' = r — — 50 x r where 
6r is an infinitesimal translation and (50 is a vector giving the orientation of the 
axis of rotation and the magnitude of an infinitesimal angle of rotation about this 
axis. We want to approximate /(r') in terms of the V/ and small translational 
and rotational variables. For r' = r + d we may always write: 

/(r') ^ /(r) + d • V/(r) (C.l) 

hence 

/(r') ^ /(r) - {6r + 50 x r) • V/(r) (C.2) 

which can be written as 

/(r') ^ fir) 6r ■ V/(r) - 50 • r x V/(r). (C.3) 

Therefore for small translations only we may write the time- varying net receive 
field as 

CUrp,tn) ^ C,,.v(rp,0) -5r(?„) • VC,,,(rp,0) (C.4) 
where we have assumed that (0) = J, the identity operator. 
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